arXiv: 1504.00013v2 [cond-mat.str-el] 20 Sep 2015 


Compressibility as a probe of quantum phase transitions 
in topological superconductors 

David Nozadze and Nandini Trivedi 
Department of Physics, The Ohio State University, 

191 W. Woodruff Avenue, Columbus, OH 43210, USA 
(Dated: September 22, 2015) 

The non-Abelian statistics of Majorana fermions, their role in topological quantum computation, 
and the possibility of realizing them in condensed matter systems, has attracted considerable at¬ 
tention. While there have been recent reports of zero energy modes in single particle tunneling 
density of states, their identity as Majorana modes has so far not been unequivocally established. 
We make predictions for the local compressibility fti oc , tuned by changing the chemical potential p 
in a semiconducting nanowire with strong spin-orbit coupling and in a Zeeman held in proximity 
to a superconductor, that has been proposed as a candidate system for observing Majorana modes. 
We show that in the center of the wire, the topological phase transition is signaled by a divergence 
of Kioc as a function of p which is an important diagnostic of the topological phase transition. We 
also find that a single strong impurity potential can lead to a local negative compressibility at the 
topological phase transition. The origin of such anomalous behavior can be traced to the formation 
of Andreev bound states close to topological phase transitions. Measurable by a scanning electron 
transistor, the compressibility includes contributions from both single particle states and collective 
modes and is therefore a complimentary probe from scanning tunneling spectroscopy. 

PACS numbers: 71.10Pm, 03.67.Lx, 74.45.-l-c, 74.90.+n 


The search for Majorana fermions in condensed mat¬ 
ter systems and their possible application for topologi¬ 
cal quantum computation 1-5 has led to several promis¬ 
ing proposals for practical realizations of Majorana 
fermions both in one (ID) and two-dimensional (2D) sys¬ 
tems. Majorana fermions can emerge in systems, such 
as topological insulator-superconductor interfaces 6,7 , 
quantum Hall states with filling factor 4 5/2, p -wave 
superconductors 8 , semiconductor heterostructures 9 ’ 10 , 
half-metallic ferronragnets 11,12 and ferromagnetic metal¬ 
lic chains 13 . As shown by Kitaev 14 , Majorana fermions 
can emerge at the ends of ID spinless p -wave supercon¬ 
ducting chain when the chemical potential is in the topo¬ 
logical regime. 

A realization of the Kitaev chain based on a quantum 
nanowire made of a semiconductor-superconductor hy¬ 
brid structure has been proposed 9,10 . In the presence of 
Rashba spin-orbit coupling, the parabolic bands for the 
two spin projections get separated. In addition, a Zee- 
man field h opens up a gap leading to an effectively spin¬ 
less ID system when the chemical potential p lies in the 
Zeeman gap. The proximity induced superconductivity 
with a gap A can result in the topological phase 9,10,15,16 . 
In this regime, the wire can be realized as a Kitaev 
chain and should have two Majorana localized zero en¬ 
ergy modes at the ends. The nanowire can undergo a 
quantum phase transition from a topologically trivial su¬ 
perconducting phase to the topological one (or vise verse) 
by changing the chemical potential or the magnetic field. 

There have been recent reports of observations of Ma¬ 
jorana fermions in tunneling and the fractional Joseph- 
son effect 16-19 . Ref. 20 has reported significant progress 
in creating the Majorana states where spatial location 
of Majoranas are detected using a scanning tunneling 


microscope. All these experimental observations of the 
existence of Majorana fermions assume that the system 
is in the topological phase and attribute the zero energy 
density of states to the proposed Majorana modes. How¬ 
ever, since there could be several other sources of zero 
bias anomaly 21-23 , the existence of the Majorana modes 
has so far not been unequivocally established. 

We propose here a definitive method to determine 
whether or not the nanowire is in the topological or triv¬ 
ial state through measurements of the gate-tuned local 
compressibility. We find very different behavior of the 
local compressibility at the edge of the wire and in the 
center. While the edge harbors Majorana modes which 
show a zero bias anomaly, the density of states in the 
center is fully gapped. The compressibility on the other 
hand shows a sharp singularity in the center of the wire 
at the topological to trivial phase transitions (see Fig. 
1). This is one of our central results. Another key result 
is that a single impurity can dramatically change the lo¬ 
cal response: a strong impurity leads to the formation of 
Andreev bound states and surprisingly this results in a 
local negative compressibility with a dip at the topolog¬ 
ical phase transitions as the chemical potential is tuned. 
An extra peak associated with the bound state appears 
in the local compressibility above the transition. 

Model and methods: We consider a ID tight-binding 
Hamiltonian 

H= 'y ' (p — Vi)<? i(J Ci a — y ' ti(c\ a Ci+la + h.c.) 

i,a icr 

+ Hso + Hz + Hint ) ( 1 ) 

where c] (c^) is the creation (destruction) operator for 
an electron on a site i , ti the nearest-neighbor hop- 
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FIG. 1. (Color online), (a) The local particle density ni oc (dashed line) and the local compressibility Ki oc (solid line) versus 
chemical potential fi/t in the center of wire. The compressibility has sharp peaks at the transition between trivial (tS) and 
topological (TS) superconducting phases. The singularity in ki oc is weakened at the edge of the wire (b). (c,d) The local density 
of states N\ oc {lj), measured relative to the chemical potential /x, in the trivial (/x = 0, blue circle) and topological (/x = 1.54, 
red star) phases. Results presented for chain length N = 256 at T = 0 with slightly broadened 5-functions in JVioc(w). 


ping and /j, is the chemical potential. The spin-orbit 
coupling and the Zeeman field terms are given by 

Hso = ^E w a i c\ +la (ia y )<j C ,'Ci' J i + h.c. and = 


_ ^S JC r c !+i 0 -( cr -)<rcr' c io-', respectively. Parameters on re¬ 
fer to the Rashba spin-orbit coupling and h to the Zee- 
man field. Vi is the on-site impurity potential. The in¬ 
teraction term H\ nt = —U ]Tb where U is the 

pairing interaction. In the clean limit a.i = a, ti = t and 

v z = o. 


We solve the model in Eq. (1) within the Bogoliubov- 
de Gennes (BdG) self-consistent approach and calculate 
the local particle density ni oc and the local density of 
states iVioc(w); (see Supplement for more details). Even 
though this is a one-dimensional problem we are justified 
in ignoring the quantum fluctuations, primarily because 
the system is proximity coupled to a bulk superconductor 
which damps out the fluctuations. 

As shown in Ref. 24, the model in Eq. (1) has 


several different phases: trivial superconducting phase 
(tS), topological superconducting phase (TS) with Ma- 
jorana fermions at ends of chain, Fulde-Ferrell-Larkin- 
Ovchinnikov phase (FFLO) with spatially oscillating or¬ 
der parameter A and non-zero magnetization, insula¬ 
tor phase (INS) with finite energy gap and normal gas 
(NG) phase without pairing and energy gap (see Fig. 2 
(b)). We are interested in a particular slice of the 
phase diagram in order to investigate the behavior of 
the compressibility across the topological to trivial phase 
transitions. We consider the Zeeman field h = 1.354 
and spin-orbit coupling a = t that shows two tran¬ 
sitions from the topological to the trivial phases as a 
function of the chemical potential /i at and /i + . 
We find the quasiparticle excitation energy E±(k) = 

el + a 2 sin 2 (fc) + ft- 2 + A 2 ±2y' / ft- 2 (e^ + A 2 ) + a 2 sin 2 (fc)e|, 
where efc = — 2fcos(fc) — /x. When h > 0, the gap 
(FA = 0) closes at |/x ± | = (2t ± Vh 2 — A 2 ), which cor- 























3 



FIG. 2. (Color online), (a) Proposal for measuring local compressibility using an SET (single electron transistor) in the 
suggested set-up 9 ' 10 of a nanowire (NW) with spin-orbit coupling in proximity to a superconductor (SC) with an applied 
magnetic field, (b) Phase diagram of ID spin-orbit coupled superconductor as function of Zeeman field h/t and the chemical 
potential p/t 24. Five different phases can be identified: trivial superconducting (tS), topological superconducting (TS), FFLO, 
normal gas (NG) and insulator phase (INS), (c) The band structure for a wire with spin-orbit coupling in a magnetic field. As 
attraction is turned on, different pairing symmetries emerge depending on the location of p (the Bogoliubov bands have not 
been shown). Within the first band, the system is described by the Kitaev model that captures the transition from the trivial 
p to topological p-wave SC. Once the second band is crossed, both interband s-wave and intraband p— wave channels become 
operative. 





FIG. 3. (Color online). Panels (a,b) The local compressibility ki oc as a function of the chemical potential p/t for link defects 
measured on either side of the link. Panel (a) Positive link defect (spin-orbit coupling ao ^ a and hopping parameter to ^ t 
different from the reference values shown for three cases (i) No defect (black dotted line); (ii) ao = to = 0.3 1 (blue thick solid 
line); (iii) cut wire with ao = to = 0 (blue thin line). The singularity in reioc at the topological phase transitions at p = p± is 
weakened in the presence of the link defect. Panel (b) Negative link defect with to = —t (blue solid line). The sharp peak in 
Kioc within the topological phase arises because of the formation of a zero-energy bound state on either side of the link defect, 
(c) Local density of states (LDOS) for a negative link defect showing a bound state at zero energy at p = 2t (blue) and away 
from zero p = 2.3t (red). 


responds to the phase transitions. The system is in the 
topological phase when p~ < \p\ < p + and in the trivial 
phase for p < —p + or p > p + (see Fig. 2 (c)). For the 
parameters chosen, p~ ss 0.84f and p + ss 3 1 . It is im¬ 
portant to note that the system is effectively a “spinless” 
p- wave superconductor so long as the chemical potential 
crosses only a single band and the transitions are from 
the topological phase to a trivial p-wave phase. Once 
both bands are crossed, the trivial phase has contribu¬ 
tions from both s (inter-band) and p- wave (intra-band) 
paring channels 5 . 

Compressibility at the topological phase transition: From 
the self consistent BdG solutions of Eq. (1), we obtain 
the local particle number n\ oc and its dependence on 
the global chemical potential p yields the compressibil¬ 
ity Ki oc = dn\ oc /dp. We find that in the center of the 


wire the local compressibility shows a logarithmic diver¬ 
gence arising from the gap closing linearly between the 
topological to the trivial phase transitions (Fig. 1 (a)). 
In contrast, the singularities are weakened on the edge 
of the wire (Fig. 1 (b)). It is useful to contrast the com¬ 
pressibility which captures the single particle and the pair 
(collective modes) density of states (DOS), from the be¬ 
havior of the single particle density of states N\ oc (oj). 
As seen in Fig. 1 (c,d), the DOS shows a gap in both the 
topological and the trivial phases except at the transition 
where the gap gets closed. However, the compressibility 
is non-zero inspite of a single-particle gap because of the 
contribution from the pairs. 

This is one of our central results. It highlights the 
fact that by measuring both the local tunneling DOS at 
the edge of the wire and the local compressibility at the 





































4 



FIG. 4. (Color online), (a) The local particle density n\ oc and Ki oc for a local potential defect Vo > P + ■ Close to the topological 
phase transitions p = /x , ki oc becomes negative. In addition, extra peaks (shown by circles) appear in the topological and 
trivial phases, (b) Density plot of Ki oc in the p/t - Vo /t plane. (c,d,e) Local density of states for various values of Vo and p. 
(c) shows how a in-gap state appears in the presence of an on-site impurity in the trivial phase (/x = 0 < p~) with s and p 
wave parings, (d) shows how a non-trivial bound state forms in the topological phase close to the topological phase transition 
at /x = 0.9£ > p~ as impurity strength Vo increases. The bound state starts to be formed when Vo > p~ . (e) The bound state 
in the topological phase for p = 1.2 1. 


center of the wire as a function of /x it is possible to un¬ 
equivocally determine when the wire is in the topological 
phase with Majorana modes localized at the edges. As a 
control, p can be varied to bring the wire into a trivial 
phase with a finite compressibility and a gapped single 
particle DOS. 

We next discuss the effect of a weak link and a single 
potential disorder on the local compressibility. 

Weak link : In the presence of a weak link, defined by a 
hopping to ^ t , the local particle number n\ oc becomes in¬ 
homogeneous. We calculate the local compressibility on 
either side of the link defect ki oc = dn\ oc /dp by differ¬ 
entiating it with respect to the global chemical potential 

p. 


For a positive link defect, i.e. to has the same sign 
as f, we find that the peaks in the local compressibility 
are weakened at the transition point (Fig. 3 (a)). In the 
limit to = 0 the wire is cut, the singularity in «i oc is com¬ 
pletely suppressed, though the compressibility remains 
finite. We also consider the case of a link with negative 
hopping parameter i.e. to = —t. Such a negative link 
can be produced by a local 7r-junction. The behavior of 
Kioc as a function of p is remarkably different from the 
positive defect. In the topological phase a sharp peak ap¬ 
pears in the local compressibility measured on the either 
side of the link (Fig. 3 (b)) at a particular value of p 0 ; for 
the chosen parameters, po = 21. This is due to the for¬ 
mation of a zero-energy bound state (Fig. 3 (c)) at po. It 
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is important to note that this zero-energy bound state is 
formed only in the topological phase. At the same time, 
it does not correspond to a Majorana mode based on the 
structure and symmetries of the corresponding eigenfunc¬ 
tions (see Supplement for more details). For p ^ p, 0 the 
bound state moves away from zero energy and no longer 
contributes to the singularity in the compressibility. 
Local potential: In the presence of an on-site impurity Vo, 
there are several interesting features in the behavior of 
the local particle density and compressibility as shown in 
the density plot in Fig. 4 (b). 

(i) Negative local compressibility: A repulsive potential 
Vo > pT can have a non-trivial effect on the local density 
and compressibility. The local particle density is found 
to decrease around the topological phase transition at 
/i = pT even as p increases. Correspondingly, the local 
compressibility ki oc becomes negative and shows a dip at 
the transition. For impurity strength somewhat larger 
than Vo > /i + , in addition to the dip at p = p~, a second 
dip appears at p = p + where also the local compress¬ 
ibility Kioc becomes negative (Fig. 4 (a)). The plot for 
the particle density shown here is schematic for clarity, 
the actual data with more details can be found in the 
supplementary information. The reason for the decrease 
of the local density and the corresponding negative local 
compressibility is tied to the formation of an impurity 
bound state (BS) above zero energy that starts to form 
close to the topological phase transitions. 

(ii) Bound states: As seen in Fig. 4 (c), In the pres¬ 
ence of an on-site impurity, the peaks at the gap edge 
are suppressed and the gap size is reduced. This can 
be understood from the fact that the trivial phase has 
both s and p-wave paring, and disorder affects the p -wave 
component more drastically than the s-wave component; 
however the spectrum remains gapped. In the topolog¬ 
ical phase where the system is effectively a “spinless” 
unconventional (p-wave) superconductor, a bound state 
is formed due to the sign change of the order parameter 
in this unconventional superconductor. 

Figure 4 (d) shows how the zero-energy bound state 
starts to form when Vo > p~. For a fixed Vo as p in¬ 
creases, the bound state becomes sharper and moves to 
zero energy. At this point the zero-energy bound state 
is detectable as an additional feature shown by a circle 


in fooc (Fig- 4 (a)). With further increase of p the BS 
moves below the chemical potential (Fig. 4 (e)). Simi¬ 
larly, a zero-energy BS forms also in the trivial p-wave 
phase for the impurity strength Vo > p + . 

For a negative impurity potential, the BS forms below 
the Fermi level and more states shift below the Fermi 
energy to enhance the local density for all p. In contrast 
to the scenario of the positive impurity potential, the BS 
does contribute to the local particle density for a negative 
impurity. As the result, the local particle density starts 
to increase as p decreases, until a sharp BS is formed. 
This once again causes the local compressibility to be¬ 
come negative around the topological phase transition; 
(see Supplement for more details). 

Conclusions: Our theoretical proposals based on the 
compressibility, in conjunction with scanning tunneling 
spectroscopy, are powerful diagnostics for detecting topo¬ 
logical phase transitions in ID spin-orbit coupled super¬ 
conductors. Specifically in the presence of local defects, 
the local compressibility can be measured using single¬ 
electron transistor (SET) spectroscopies 25 . Ref. 25 has 
in fact used the SET in a different context to measure 
the inverse compressibility locally on a graphene sam¬ 
ple as a function of the back-gate voltage or carrier den¬ 
sity. We expect the same technique can be applied to the 
spin-orbit coupled nanowires - superconductor devices to 
detect the topological phase transition guided by our pre¬ 
dictions. 

Some of the most promising directions to experimen¬ 
tally investigate are: (a) the sharp peak in the compress¬ 
ibility at the topological phase transition tuned by the 
Zeeman field in the clean wire, and, (b) the negative 
compressibility induced by the on-site impurity in the 
topological phase. In general it will be useful to see the 
interplay between local scanning and local compressibil¬ 
ity spectroscopies for giving insights into single particle 
and collective modes. 
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In this supplement, we present: (I) details of the 
Bogoliubov-de Gennes (BdG) mean field theory for the 
nanowire in proximity to a superconductor; (II) informa¬ 
tion about the eigenstates associated with zero-energy 
bound states that are localized around a negative link 
(7r phase shifted) in the topological phase but are not 
Majorana modes; (III) additional information about the 
local particle density, local compressibility and the local 
density of states (LDOS) in the presence of an on-site 
impurity. 

I. Bogoliubov-de Gennes (BdG) approach: 

Starting with the full Hamiltonian in Eq. (1) in the 
main text, we perform a mean-field decomposition of the 
interaction term in the s-wave local pairing amplitude 
channel, A; = —U (cj^cj^) . In addition, we can add an 
on-site impurity potential or a weak link. The mean field 
quadratic Hamiltonian is given by: 

H = ~ - V i) C h C i° - X ^(4xG+1<t + h.C.) 

i,cr io 

+ ??so + Hz + Ai(cj|c|^ + h.c.). (1) 

i 

We next diagonalize Eq. (1) using the Bogoliubov 
transformation c ia = J2 n ( TnU nlJ {i ) - a^ n v na {i)) which 
leads to BdG equations: 



a 

0 

A ^ 

( Un t(j) 

—a 


-A 

0 

Unl(j) 

0 

-A* 

~H t 

—a 

1 -Vntti) 

A* 

0 

a 

-Hj 

\ v n ±(j) 


/ «nt(j) 
Unl(j) 

I -Wt(j) 

V v ni {j) 

( 2 ) 


where the excitation eigenvalues E n > 0. Here 

Hfj'lLfiuij'') — ( E fl ~\~ Vi^Uncr (^) ^(^ncr(^ “hi) - ! - 

Una{i - 1)), au nt {i) = a(u ni (i - 1) - u ni (i + 1)) and 
A u na {i) = AiU na (i). In terms of the eigenfunctions of 
the BdG hamiltonian, the local pairing amplitudes and 
local densities are given by: 


Ai = -U y ^ u n ^(i)v n ^(j) tanh(E' n /2T). (3) 

En> 0 


and 


n i = \E [\ u ™(i)\ 2 .f( E n) + K ct (*)| 2 (1 - f{E n ))] , 

n,a 

( 4 ) 


where f{E n ) is the Fermi function. The solutions of 
the BdG equations are iterated until self consistency is 
achieved at each site. In these studies we have set T = 0. 
In the problem at hand, since the pairing in the nanowire 
is imposed by proximity with a bulk superconductor, it 
may be not necessary to obtain the pairing amplitudes 
self consistently, however in the presence of disorder (site 
or link), the self consistency becomes more relevant and 
is included in our investigation. 

Once self consistency is achieved the single-particle 
density of states is obtained from 

Ni(u) = ^2 [\una{i)\ 2 &(u - E n ) + \v na ((i)\ 2 S(u) + E n )] . 

n,cr 

( 5 ) 

In the text the local density from Eq. (4) is denoted 
by ni oc and the local single-particle density of states by 

N\ oc(w). 

II. Zero-energy bound state due to a weak link: 

As shown in the main text a negative link impurity to = 
—t induces zero-energy bound states in the topological 
phase. Are these Majorana modes or simply Andreev 
bound states? In order to get some insight about these 
inodes we look at the corresponding eigenfunctions u a (i) 
and v a (i). As seen in the Fig. 1 (a), close to link impurity, 
the eigenstates do not satisfy the symmetry conditions 
u a {i) = v*(i) or u a {i) = that is required for a 

Majorana mode (Fig. 1 (b)). Thus, even though these 
two zero-modes are in the topological phase these are 
not Majorana modes. 

III. Local particle density and compressibility: 

We present here the actual data in Fig. 2 of the supple¬ 
ment for the local particle density and the local com¬ 
pressibility in the presence of an on-site impurity po¬ 
tential ji + > Vo > n~, corresponding to Fig. 4 (a) 
of the main text. Close the topological phase transi¬ 
tion at /.t = /W = 0.84t, the local particle density starts 
to decrease even as the chemical potential increases and 
the local compressibility ki oc becomes negative ( Fig. 2 
(a,b)). Also, an extra peak appears in k\ oc in the topo¬ 
logical phase (Fig. 2 (b)) due to the formation of an An¬ 
dreev bound state as explained in main text. For an 
on-site impurity potential Vo > U + , close the topological 
phase transitions yu. = /x =l= = 3f, the local particle den¬ 
sity once again starts to decrease as the chemical poten¬ 
tial increases and local compressibility becomes negative 
(Fig. 2 (c,d)). The additional extra peak in the com¬ 
pressibility now appears in the trivial phase in addition 
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FIG. 1. (Color online), (a) The bound state eigenfunctions at zero-energy in the presence of a negative link impurity to = —t 
between sites 15 and 16. (b) The bound state eigenfunctions at zero-energy for a broken link to = a = 0 between sites 15 an 
16. From the symmetry conditions u CT (i) = — v a (i) required for a Majorana fermion, it is evident that while the bound state in 
(b) is a Majorana bound state that is topologically protected to be at zero energy, the bound state in (a) is an Andreev bound 
state that happens to be at zero energy. The system size N = 256, and parameters a = t, h = 1.35t at T = 0. 
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FIG. 2. (Color online). The local particle density n\ oc and the compressibility ki oc versus chemical potential /r/f for two values 
of the on-site impurity potential Vo- 


to peak in the topological phase (Fig. 2 (d)), once again 
arising because of an Andreev bound state. 

As discussed in the main text, we discussed the role 
of an Andreev the bound in the presence of an on-site 
impurity potential and how that results in a decrease of 
the particle density with increasing /r. This was shown 
previously in the LDOS for a positive potential Vq > 


0. Here, we present the results for the LDOS in the 
presence of a negative potential V 0 < 0. In the latter 
case, the local particle density starts to increase as the 
/1 increases because the bound state now forms below the 
Fermi level so more states shift below the Fermi level 
(Fig. 3). Thus in contrast to the scenario of the positive 
impurity potential, the bound state does contribute to 
the local particle density for a negative on-site impurity. 
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FIG. 3. (Color online). Formation of the bound state close 
to the topological phase transition for a negative on-site po¬ 
tential. 










